Transport equation: extension of classical results for div b 6 BMO 
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Abstract. We investigate the transport equation: u t + b-\7u = 0. Our result improves 
the criteria on uniqueness of weak solutions, replacing the classical condition: div b G 
3 ■ by div b £ BMO. 
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Oh ■ 1 Introduction 

< 

The goal of this note is to improve classical results concerning the Cauchy problem 
for the transport equation. The basis of our analysis is the following system 

u t + b-Vu = in M n x(0,T), 
u\ t =o = u on R n , 

> . 

where u is an unknown scalar function, b - some given vector field and Mo is an initial 
datum. 

The transport equation is one of the most fundamental examples in the theory of par- 
tial differential equations. It describes the motion of matter under influence of the velocity 
field b. Classically, for smooth data b and Uq, (1.1) is solvable elementary by the method 
of characteristics. In the language of the fluid mechanics, (1.1) says that u is constant 
along streamlines defined by the Lagrangian coordinates. This physical interpretation 
gives enough reasons for (1.1) to be intensively studied from the mathematical point of 



O 

as 



view. Here we want to concentrate on the optimal/critical regularity of the vector field b 
to control the existence and uniqueness of the solutions. 

The classical results require that the vector field b must fulfill 

(1.2) divfoe J L 1 (0,T,L oo (R")) ) 

then we are able to obtain existence and uniqueness of weak solutions to (1.1). The 
objective of our investigations is to relax the condition (1.2) to 

(1.3) div 6 G Lx(0,T; BMO(M. n )). 

Let us observe that this "slightly" broader class than (1.2) is on the boundary of known 
counterexamples [6]. For any p < oo we are able to construct such b e VPp(lR n ) (time 
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independent) to obtain an example of the loss of uniqueness to (1.1). On the other hand 
the L?MO-space appears naturally in many considerations, since it is the limit space 
for the embedding W^(R") C B MO (R n ), where the Loo-space is not reached. We are 
able to prove existence and uniqueness of weak solutions to (1.1) in the case of bounded 
solutions and improve the uniqueness criteria for L p -solutions. Additionally we show 
a result concerning stability with respect to initial data. Our approach follows from 
techniques introduced in [11] to improve the uniqueness criteria for the Euler system in 
bounded domains. 

Results fundamental for our issue have been stated by R.J. DiPerna and P.L. Lions in 
[6], where general questions concerning the well posedness of the problem found positive 
answer under condition (1.2). An interesting extension of the theory has been made by 
L. Amrozio [1], for the case of bounded solutions replacing the condition b G W±(R n ) 
by b G BV(M. n ). In the literature one can find also numerous works on generalizations 
of the mentioned results on broader class of function spaces [2] , [4] , [7] , [8] , [9] , but positive 
answers still require condition (1.2). 

In the present note we consider weak solutions meant in the following sense: 

We say that u G Ax>(0, T; L p (R ra )) is a weak solution to (1.1) iff the following integral 
identity holds 



(1.4) / / u(f) t dxdt+ / / div b u<pdx + / / b-V(j)udxdt = — I u (f)(-,0)dx 
Jo Jk™ Jo Jw n Jo Jm.™ it" 

for each G C°°([0, T];C X (W 1 )) such that <p\ t=T = 0. 

Let us state the main results of this paper. First we start with the case of pointwise 
bounded solutions, in that case our technique delivers the most complete result. 

Theorem A. Let T > 0, b G Li(0, T; W^ oc) (R™)) ; u G L 00 (R™) ; additionally we 



assume 



(1.5) div b G Li(0, T; BMOiW 1 )), & G Li(0, T; Li(R")) 



(1.6) and supp div b(-, t) C B(0, R) for a fixed R > , 

where B(0,R) denotes the ball centered at the origin with radius R. 

Then there exists a unique weak solution to the system (1.1) such that 

(1.7) u G L oo (0,T;L oo (R n )) and ||u||l O o(0 ) T;L oo (r»)) < ||uo|Uoo(r» 



The above result guarantees not only the uniqueness of solutions, but also their exis- 
tence. It is a consequence of a maximum principle, which is valid for the Loo-solutions. 
To show (1.7) condition (1.5) is not needed. The main difference to the classical results 
[6] is that having (1.2) we are able to construct the L p -estimates of the solutions for finite 
p. In our case the condition (1.3) is too weak to get such information. Additionally we 
are required to add an extra condition (1.6), which is the price of our improvement of this 
classical criteria. 
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The next result concerns stability of solutions obtained in Theorem A with respect to 
perturbations of initial data in lower spaces. 

Theorem B. Let 1 < p < oo and b fulfill assumptions of Theorem A. Let u ,Uq e 
Loo(lR n ) and (uq — u ) e L p (R n ) such that sup fcgN ||-UolUoo(K n ) + || u o|Uoo(K") = m < 00 an< ^ 
(uq — u ) — > in L p (R a ) as k — > 00. Then 

(1.8) (u k -u)^0 inL oo (0,T;Lp(R n )) as fc -> 00. 



The last result concerns the uniqueness criteria for L p -so hit ions to (1.1). 

Theorem C. Let 1 < p < 00, - + ± = 1, b e Li(0, T; W^ oc) (R n )) and conditions 
(1-5), (1.6) be fulfilled. Let -a 1 , -a 2 &e two weaA; solutions to (1.1) with the same initial 
datum and u 1 , u 2 e L^O, T; L p (R n )); then u 1 = u 2 . 

The above three results are proved by a reduction of considerations to an ordinary 
differential inequality of the form 

(1.9) x = x\ax, x\t=Q — 0. 

The Osgood lemma yields the uniqueness to the system (1.9). This observation forms 
our chain of estimations to have a possibility to adapt information obtained from the 
Gronwall inequality. Due to low regularity of solutions, our analysis requires a special 
approach. The main tool, which enables us to show the main inequality in the form of 
(1.9), is Theorem D stated below. 

Theorem D. Let f e BMO{R n ), the support off be bounded in R n and g e Li(R n ) n 
L oc (R n ), then 

(1-10) I / fgdx\ < C ||/||bmo(r»)IMUi(r») [I ln IMUi(R")l + ln ( e + llfi , IU o(K™))] • 

The above inequality can be viewed as a representant of the family of logarithmic 
Sobolev inequalities [3], [5], [10], however there is one important difference between this one 
and others, here an extra information about derivatives of the function is not required, in 
contrast to L^ — BMO inequalities. The crucial assumption is the boundedness of the 
support of the function /, it is a consequence of results of the classical theory [13], [14] 
and unfortunately it is not expected that it could be possible to omit this restriction. 
Methods of proving (1.10) distinguish this result from others, too. They base on relations 
between the Zygmund space L ln L and Riesz operators. Theorem D has been proved 
in [11], applied successfully to the evolutionary Euler system. Outlines of the proof of 
Theorem D one can find in the Appendix. 

The below remark shows us a possible generalization of stated theorems. 

Remark. The results stated in Theorems A, B and C can be easily extended on the 
following linear system 

, . u t + b • Vw = cu + / in Qx(0,T), 

u\ t= o = u on f2 
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in an arbitrary domain Q C M n with a sufficiently smooth boundary dVL, say C 0+1 , enough 
to allow integration by parts, and with given 

c,f G ^(O^L^R")) and b-n = on dQ x (0,T), 

where n is the normal vector to the boundary dfl. 

Additionally, we find a natural generalization of (1.5)- (1.6) 



div b = + H k such that 

k=l 

(1.12) G L^O, T; Loo(fi)), H k G Li(0,T; BMO{Q)) and 

oo 

y2 II^IUi(o,T;BMO(Q)) < oo with sup diam supp H k < oo. 

fc=l fceN 

In the case of bounded Q condition (1.12) simplifies itself and (1.6) is automatically 
satisfied. We leave the proof of Remark to a kind reader, it is almost the same as for 
(1.1), the estimations are just more technical, but the core of the problem is the same. 

Thoughout the paper we use standard notation. L p (R n ) denotes the common Lebesgue 
space, generic constant are denoted by C. Let us recall only the definition of the BMO- 
space. We say that / G BMOiW 1 ), if / is locally integrable and the corresponding 
semi-norm 

(1-13) ||/||bmo(m») = sup / \f(y) - {f} B ( x ,r)\dy 

xm n ,r>0JB(x,r) 

is finite, where {f}B(x,r) = \B(l,r)\ $B(x,r) f(v) d y and B { x , r ) is a ball with radius r centered 
at x - see [12]. The above definition implies that (1.13) is a semi-norm only, however in 
our case from assumptions on div b follows 

II div6|| Ll(K n) < \supp div b\\\ div b\\ B MO(R™) 
which is a consequence of the properties of the support restricted by (1.6). 

2 Proof of Theorem A 

Our first aim is to prove the existence of weak solutions to (1.1). To construct them 
we find the following sequence of approximation of the function b and initial datum uq. 
We require that ¥ G C 00 ^™ x (0,T)) and supp div ¥(■, t) C B(0,2R) and b e -> 
b in Li(0,T; Wh loc JM. n )) with suitable behavior of the norms. For the initial datum we 
find Uq G C£°{R n ) with u e ^* u in L 00 (R") as e -> and 

(2- 1 ) ll M olUoo(K") < HuolUooCK")- 

Then we consider the following equation with smooth coefficients b e and initial data u e . 

, , u\ + b e ■ Vu e = in R n x (0, T), 

^ • ' u% =0 = Un on R n . 
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The method of characteristic implies the existence of smooth solutions to (2.2) for 
t G (0, T). Omitting the characteristic coordinates system we find an e-independent 
estimate for the solutions by a simple Mouser's technique. For any even p we easily get 

-— / u ep dx-- [ dlvb e u ep dx = 0, 



pdt J p 
which implies 

II €11 / II en - fn lldiv b e \\r (m n)ds 

Since div u e is smooth and e > is fixed, hence letting p — > oo, we get immediately 

( 2 - 3 ) 1 1 ^ £ 1 1 Loo (0,T;Loo (M™)) < |KIUoo(R")- 

Note that we do not use any uniform bound on div6 e . 

Now we pass to the limit with e — > in (2.2). The solutions to (2.2) are classical, in 
particular it implies they fulfill the following integral identity 

(2.4) -f f u £ (j) t dxdt- I [ div b e u £ (f)dx- I [ b e -V^u'dxdt = [ u £ <f)(- , 0)dx 

Jo JR n Jo Jr" Jo Jr 11 Jr™ 

for any G C°°([0,T]; C~(R n )) such that (p\ t=T = 0. 

Estimates (2.1) and (2.3) imply that for a subsequence e& — > 

(2.5) w efe ^* « in Loo(0, T; (M n )) with ||w||l 00 (o,t ; l 00 (r»)) < IKI|loo(R»))- 

Then taking the limit of (2.2) for — >• 0, by properties of sequences {b e } and {u € }, we 
obtain 

(2.6) — / / u(f) t dxdt — / / div bu<f)dx— / / b-V<pudxdt= I u (f)(-,0)dx 

Jo Jr u Jo Jr u Jo Jr u Jr u 

for the same set of test functions as in (2.4). To simplify the notation we will write e — > 
instead of — > 0. 

This way the limit m, defined by (2.5), is a weak solution to (1.1). It is clear that 

(2.6) it is its distributional version - (1.4). However the high regularity of test functions 
required in (2.6) does not allow us to obtain any information concerning the uniqueness 
of solutions to (2.6) in a direct way. To solve this issue we start with an application of 
the standard procedure. We introduce 

(2.7) S £ (f) = m £ *f= [ m £ (--y)f(y)dy, 

Jw n 

where m e is a smooth function with suitable properties tending weakly to the Dirac delta 
- see (5.7) in the Appendix. Applying the above operator to (1.1) we get 

(2.8) d t S e {u) + S e (b ■ Vu) = 0, 

where b ■ Vu = div (bu) — u div b and the r.h.s. is well defined as a distribution. In fact 
(2.8) implies that d t S e (u) is well defined as a Lebesgue function, too. 
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\x\ 


< 1 




(2.10) 




[4 l] 


1 < 


\x\ < 2 


with V7r r < 










> 2 





We state equation (2.8) as follows 

(2.9) d t S e (u) + b-VS e (u) = R e , where R t = b ■ VS e (u) - S e (b ■ Vu). 

Standard facts follow (see (5.8) in Appendix) the remeider is controlled in the limit: 
R t — > in Li(0, T; L 1( 7 oc )(IR n )). Since i? e convergences locally in space, only, we introduce 
a smooth function n r : W 1 — > [0, 1] such that 7r r (:r) = 7Ti(^) and 

r 

In order to prove the uniqueness for our system it is enough to consider (2.9) with zero 
initial data (due to its linearity). Since we are forced to localize the problem, we multiply 
(2.9) by S e (u)n r and integrate over the space, getting 

-— / (S e (u)) 2 7T r dx — - [ div b (S e (u)) 2 'n r dx — — f b ■ Vn r (S e (u)) 2 dx 
(2 11) J Rn •' IRn 

= / R e S e (u)ir r dx. 
Jm™ 

Then integrating over time, using properties of S € and letting e — > 0, next differentiating 
with respect t we obtain 

X d /* X /* X /* 

(2.12) -— / u 2 T\ r dx — - I div bu 2 ir r dx — ^ I b ■ VK r u 2 dx for r > 0. 

^ Jr™ ^ it" ^ </R" 

The r.h.s. of (2.12) is estimated as follows 
(2.13) 



/ b- Vn r u 2 dx\ < C\\u\\ 2 Loo( ^ n) 





b\ 




1 + 


\x\ 



1 + 1x1)1 Vn r \dx — > as r — > 00. 



0. 



*\B(0,r) 

By definition (1 + |x|)|V7iy| < C, because the support of V7r fc is a subset of the set: 
{r < |x| < 2r}. By (2.5) the norm IHIl^ is controlled, too. Then letting r — > 00 in 

(2.13) , we get 

~\. d f X /* 

(2.14) -— / u 2 dx — - I divbu 2 dx = with w|t=o 

^ Jr™ ^ JR™ 

On the other hand, the function u can be viewed as a difference of two solutions to 
(1.1) with the same initial datum, hence they satisfy the bound from (2.5), so there exists 
m > such that 

(2-15) IMUooCo.TjLooCR-)) < m. 

The application of Theorem D to the second integral in (2.14) leads us to the following 
inequality 

(2.16) 4lk 2 |Ui(iR") < c \\ div u\\bmo(R")\\u 2 \\ Li (R") [\ In ||m 2 ||li(r™)| + hi(e + m)] 

with the initial datum ||w 2 |i=o||z,i = and m from (2.15). The Osgood lemma applied to 
(2.16) yields the uniqueness to (1.1). Note that (2.16) has the form of (1.9) mentioned in 
Introduction. Theorem A is proved. 
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3 Proof of Theorem B 

The next result concerns the stability of solutions from Theorem A. We start with the 
mollified equation (2.8) for reasons same as previously, testing it now by \S e (u)\ p ~ 2 S e (u)ir r 
with p as in Theorem B. Repeating the considerations for (2.9)-(2.14) we deduce 



(3.1) \u k -u\ p dx<l | div b\ \u k - u\ p dx. 

dt J~R n Jw 1 

For a given 1 > e > 0, we fix K e G N such that for all k > K e 
(3-2) H-u \\l p < e. 

Let X p = \u k — u\ p , then by Theorem D (m as in (2.15)) (3.1) reads 

f Xdx < f | div b\Xdx 

(3.3) dt Jjgn J K n 

< C || div 6||bmo(r»)||-X'||li(r») [I In ||-X"||li(r»)| + ln(e + 2m)] , 



with J Rn X(x, 0)dx < e. 

By assumptions the r.h.s of (3.3) is at least locally integrable, hence f Rn X(x, t)dx is 
uniformly continuous. There exits a positive time T so small that 

(3.4) sup / X(x,t)dx < e" 1 . 

It follows that the function w\ \nw\ will be considered as increasing, since J Rn X(x, -)dx 
on the chosen time interval takes the values only from the interval [0, e -1 ]. Monotonicity 
allows us to introduce a function B : [0, T ] — > [0, oo) such that 

(3.5) ^S = C7 ||div6|| B MO(R»)5[|liiS|+lii(e + 2m)] and B\ t=0 = e. 

The definition of B guarantees that it is an increasing and continuous function, thus there 
exists Ti > such that < T\ < T and 

(3.6) B{t) < e- 1 < 1 for t G [0,71]. 
Taking the difference between (3.5) and (3.3) we get 

[B - Xdx] > C || div b\\ B MO(R™y 

(3.7) ■ uv r J * n . r 
■ B\\nB\- Xdx\\n Xdx\ + ln(e + 2m) {B - / Xdx) 



d 

di [ 



with B{0) - j Rn X(x, 0)dx > 0. 

Since the monotonicity of the function w\ \nw\ on [0, e -1 ] implies 

(3.8) {B\\nB\- I Xdx\\n I Xdx\){B - I Xdx)>0, 
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remembering that we consider t G [0,Ti], from (3.7) we get 

(3.9) 0< I X(x,t)dx <B(t) for te[0,7\]. 

The above fact reduces our analysis to the considerations of the function B. Additionally, 
by the choice of the time interval it follows that B(t) < 1 for t G [0, Ti], hence we can use 
the estimate 

(3.10) lln^l^lne- 1 for f e [0, 7i] . 
Solving (3.5) we get 

(3.11) B(t) <eexp|c , [ln(e + 2m) + lne- 1 ]^ /(s)dsj < Cee- c °^ fWd8 , 

where f(t) = \\ <tivb(-,t)\\BMO(R n ) an d C depends on data given in Theorems A and B. 

Next, we choose T 2 so small that < T 2 < Ti and C j^ 2 f(s)ds < 1/2, then (3.11) 
yields 

(3.12) sup B{t) < Ce 1/2 . 

te[o,T 2 ] 

Here we shall emphasize that T 2 is independent from the smallness of e - see (3.2). Thus 
we are able to start our analysis over from the very beginning, but for the initial time 
t — T 2 . Since C in (3.11) is an absolute constant we find the next interval [T 2 ,T 3 ], where 
we obtain 

(3.13) sup \\u k - u\\ Lp{R n) < Ce 1 / 4 

te[T 2 ,T 3 ] 

for all k > K e - see (3.2). Since T is fixed and finite and by the assumptions / G £i(0, T), 
we are always able to cover the whole interval [0, T] in finite steps, so finally we obtain 

(3.14) sup B(t) < Ce a . 

te[o,T] 

with a > defined by the properties of / and again C depending on all data, but 
independent from e. Letting e — > we prove (1.8). Theorem B is proved. 

4 Proof of Theorem C 

Our last result describes the uniqueness criteria for weak solutions, provided their 
existence in the /^(O, T; L p (R n )) -class. The problem reduces to (1.1) with zero initial 
data and u G Loo(0,T; L p (IR n )). To work in optimal regularity of coefficients we consider 
(2.9) 

d t S e (u) + b-VS € (u) = R e ^0 in 1^(0, T;L 1(k)c) (R")). 



8 



Next, we introduce the renormalized solution for (1.1) - we refere here to [6] where this 
approach has been developed. Take (3 G C 1 (M), i.e. H/SHl^r) + II/^'Hl^r) < oo, then 

(4.1) d t p(S e (u)) + b ■ Vf3(S e (u)) = R e P'(S e (u)) 
which implies the limit for e — > 

(4.2) dtP(u) + b • V/3(u) = 0. 
As the function j3 we choose T m : R — > [0, m p ] such that 

t a o\ /n / \ f l S ! P f° r I S I < m 

4 -3 T m (s) = i 1 L , L 

L m p tor |s| > m 

defined for fixed m G R+. T m is not a C^-function, but a simple approximation procedure 
will lead us to (4.2) with (3 = T m . 

Since we do not control integrability of all terms in (4.2) we use the function 7i r from 
(2.10) to localize the problem getting 

(4.4) — / T m {u)T[ r dx < [ | div 6| T m {u)it r dx + / \b • Vn r \T m (u)dx. 

dt Jl« jRri Jr« 

For fixed m and r letting to infinity the last term vanishes, so we obtain 

(4.5) ^ f T m (u)dx< f | div b\T m (u)dx. 

Theorem D applied to the r.h.s. of (4.5) yields 
(4-6) 

^||^m(w)||Li(R») < C|| div6|| S A^o(iR«)||r r7l (ti)|| £l ( MTl ) [| In ||T ro (u) Hl^r") I + ln(e + m p )] 

with ||T m (u)(-,0)|| £l (Rn) = 0. 

The same way as in the proof of Theorem A, the Osgood lemma yields T m (u) = 0. 
Letting m — > oo, by (4.3) we conclude u = 0. Thus, ui = u 2 - Theorem C is proved. 



5 Appendix 

Sketch of the proof of Theorem D. We proceed almost as in [11]. The assumption of 
the boundedness of supp / allows us to consider the studied integral (the l.h.s. of (1.10)) 
on a torus T n = W l /(dX a )(= [0,d) n ) with sufficiently large d guaranteeing that supp / 
can be treated as a subset of T n . Consider the Hardy space on T n with the following 
norm 

n 

(5-1) IbllwHT") = IMUi(H>) + ll-ftfc#IUi(T")> 

k=l 

where R k are the Riesz operators - [12], [13]. Since BMO{T n ) = (^(T"))*, we get 
(5-2) | / fgdx\ < ||/||BAfO(T»)IMIwi(T»)- 
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Hence to control the norm (5.1) an estimate of H-Ra^Hl^t 11 ) is required. The classical 
Zygmund's result [14] says: 

(5.3) 11-Rfc^lUi cr») < C + C J \h\\n + \h\dx, 

where ln + a = max{lna, 0} and constants C depends on d, so on the diameter of supp /. 
Let us observe that ln + (g/A) = In g — In A for g > A and 

(5.4) \\ng\ g > x \ < ln(l + ||g|| Loo ( T n)) + l In 9 | ff > A | < 2m(l + ||c/|| Loo(T ™ ) ) + | In A|. 

1 + IMIwt™) 

Taking h = — in (5.3), employing (5.4), we conclude 

(5.5) ||i2fc^|Ui(T") < C||»IUi(t») + C / \g\ [ln(l + \\g\\ Loo (T")) + I In ||ff|Ui(T»»)|] dx 

Inequalities (5.2), (5.5) yields (1.10). 

The commutator estimate. Let us recall the well known facts concerning the mollifi- 
cation of the equation and the behavior of the commutators. 
Introduce mi : W n — > [0, oo) such that 



mi (x) = N n 



exp{-j^} for \x\ < 1 
for \x\ > 1 



where the number N n is determined by the constraint f Rn m\dx — 1. Then for given e > 
we define 

. . 1 .x, r 

(5.6) m e (x) := — m\{— ) with / m e dx = 1. 

It is clear that m e — >■ 5 in D'(R n ), where <5 is the Dirac mass located at the origin of W n . 
The function m e introduces an operator S e : Lm oc ^(W n ) — > C°°(lR n ) 

(5.7) S e (h) = m e * h = / m e (x — y)h(y)dy. 

The standard theory, see [6], guarantees the following estimate for the commutator 

(5.8) b ■ VS € (u) - S e (b ■ Vu) -> in L^O, T; L 1(ioc) (M n )), 

provided that u G ^(0, T; L p(M (IT)), 6 G L^O, T; Wj (Joc) (R n )) and 1 = i + ^ for 
p G [1, oo]. The proof of (5.8) belongs to the by now classical theory, we omit it here and 
refere again to [6]. 
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